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ABSTRACT 



> 

o 

o\ . 

— ■ Context. With the superb angular resolution of the Chandra Observatory, it is now possible to detect X-ray point sources, either embedded in 
' galaxy clusters or along the cluster line of sight, which could not be resolved by previous instruments. This now allows studies of source counts 
, in distant cluster fields. 

Aims. We want to analyze the inner region of clusters of galaxies to check for the presence of any overdensity of X-ray point sources embedded 
in the gas diffuse emission. These point sources are possible AGN belonging to the clusters and could contaminate the cluster emission. 
Methods. We used a sample of 1 8 distant (0.25 < z < 1 .0 1 ) galaxy clusters from the Chandra archive to construct the log N— log S , in both the 
O j 1 soft and hard energy bands, for the X-ray point sources detected in the central cluster region to be compared with the counts of point sources 
detected in similarly deep fields without clusters. 

Results. We find a ~ 2 a excess of cluster region sources at the bright end of the log N- log 5 . The radial distribution of the brightest X-ray 
point sources confirms this excess and indicates that it is confined to the inner 0.5 Mpc of the cluster region. 
(— j Conclusions. The results suggest the possible existence of X-ray sources belonging to the cluster (most probably AGN, given their 0.5-10 keV 
• • ■ luminosity ranging from 10 43 to 10 44 erg s ): on average one every three clusters. Unlike previous studies, which have mainly investigated the 
, ^ point-source population in the vicinity of the galaxy clusters, the present study analyzes the content of point sources within the 1 Mpc region 
, covered by the cluster extent. Our work confirms the findings of other investigators who analyzed the central 1 Mpc region of more massive 
» I ■ clusters and/or groups in a similar redshift range. The X-ray source excess found here is much smaller than the excess of radio galaxies found 
C\3 recently in high-z X-ray selected clusters, possibly due to the better sensitivity of the radio observations. 

Key words, galaxies: clusters: general - galaxies: high redshift - galaxies: evolution - galaxies: active - X-rays: general - X-rays: galaxies: 
clusters 



1. Introduction 



Clusters of galaxies harbor a wide diversity of galaxy populations, so they are id eal laboratories for studying ga laxy formation 
and evolution. While there is clear evidence of the evolution in cluster galaxies dButcher & Oemlet 1978i 1984 ), the evolution 
with redshift of cluster active galactic nuclei (AGN), as well as the possible prevalence of AGN in cluster environment with 
respect to the field, is still an unresolved issue. A related open debate is whether the cluster environment plays a role in the 
probability of galaxies to develop star-forming or AGN activity. 

In the radio domain, searches for active galaxies in nearby galaxy clusters have been carried out for a long time (see among 
Qthers lOwenll975l : lFantill984T:lLedlow & Owenll995alfbT) . The Radio Luminosity Function (RLF) of nea rby cluster radio galaxie s 
was found to be statistically indistinguishable from that of the fiel d, both in shape and n ormalization dLedlow & Owen 1996 ). 
Recently, however, a radio survey of high redshift galaxy clusters (IBranchesi et al.ll2006h has provided evidence of changes in 
the RLF of the distant cluster radio galaxies as compared to the local one. All these studies show that radio galaxies in clusters 
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are strongly centrally concentrated and that their radial distribution essentially follows that of early type galaxies. As a result 
of the shape of the RLF and of the high sensitivity of radio telescopes, the typical radio luminosities found are not very high. 
Therefore the nuclear activity of these radio galaxies is not prominent at optical and/or X-ray wavelengths. On the other hand, 
studies from optical surveys seem to suggest that AGN are relatively rare in cluster environments. Cluster members that show 
evidence of AGN activity in their o ptical spectra are only 1% of all cluster galaxies, while AGN are more common (5%) in the 
field population (Dressie r et al. 1985). Furthermore, there is no evidence of any increase with redshift of the AGN fraction in 
clusters up to z ~ 5 dDressler et alJll999B . 

Since the pioneering work of He nry & Briel ( 1991), who used RO SAT data to first reveal the presence of significantly more 
X-ray point sources around Abell 2256 (z=0.06) than one would expect by chance, several studies have appeared in the literature 
that indicate X-ray point-source excesses in the vicinity of low- and high-z clusters with respect to the field. Many of these 
sources have been confirmed to be cluster members in several ca ses. At low redshift the excess seems to b e largely due to the 
low-luminosity AGN (~10 41 erg s -1 ) associated with the cluster dLazzati et al.lll998t ISun & Murravll2002l) . At higher redshift, 
more luminous (~10 42 43 erg s _1 ) AGN are found in the fields of several clusters. It is evident from these studies that there is a 
population of obscured, or at least optically unremarkable, AGN in galaxy clusters. An AGN identification may not be obvious a t 
optical wavelengths wh ere the dusty AGN o r those with weak emission lines may be not recognized as such dMartini et al ]|2002|) . 
A more recent work bv lMartini et al.l (120061) . who completed a spectroscopic survey of X-ray point sources in eight low-redshift 
clusters (z<0.3), find s that cluster galaxies host AGN more frequently than previously thought, a factor 5 higher than found by 



Dressier et al.l (11985). Thus the optical spectroscopic surveys alone can underestimate the number of AGN in clusters. This is one 



of the reasons for the recently increased popularity of X-ray and radio wavebands for AGN detection at high redshift. 

The capabilities of the current genera tion of X-ray observatories like Chandra ( van Speybroeck et alJl997 ) or XMM-Newton 
(Striider et al. 2001; Turner et al. 2001) have triggered and extended to higher r edshift these type of analyses (among others 
Cappi et al.l200ll:|Pentericci et al.l2002tlMolnar et al.l2 002: Ca ppelluti et al.l20"05 ) thereby enabling studies with very fine spatial 
details. The superb angular resolution of the Chandra observatory, as well as the high sensitivity over the full X-ray band, 
now allow the detection of those X-ray point sources within the cluster region, which could not be resolved from the diffuse 
emission with previous instruments. The works mentioned above have shown an abundance of point sources in the direction 
of individual high-z clusters as compared to the field. Different hypotheses for the apparent overdensities are examined by 
the different authors: 1) statistical variance of cosmic background sources, 2) gravitationally lensed background sources, 3) 
AGN/quasars and/or powerful starburst galaxies associated with the clusters. The last hypothesis is now considered the most 
plausible one. 



Cappel luti et al.l (120051) performed the first systematic search for X-ray point sources at the outskirts of 10 Chandra high-z 



clusters (0.24 < z < 1.2) and found a factor 2 overdensity > 2<x significance level in 40% of clusters fields. They speculate 
that the most likely astrophysical interpretation of the overdensity is that the X-ray sources are AGN that trace the filaments 
connected to the clusters. This speculation is supported by the results of lCappelluti et alj (120071) . who show that AGN and clusters 
are strongly correlated on scales betwe e n 2.5 and 10 Mpc (smaller separations were not investigated). Further support for this 
interpretation comes from lD'Elia e t al. (2004), who studied the field surrounding the cluster 3C295, which clearly exhibits a 
strong and asymmetric clustering of X-ray sources on scales of a few arcminutes. 

Most of the above-mentioned studies analyzed the whole field around the clusters. In recent years it has become evident the 
importance of the location of AGN with respect to the cluster center, which could constrain different scenarios of AGN triggering. 
The study by Ijohnson et al.l ((2003) of MS 1054-03 at z = 0.83, for instance, indicates that the X-ray AGN excess avoids the 
central regions of the cluster. Interestingly, AGN are not distributed randomly within the cluster but tend to populate the outer 
1-2 M pc, suggesting that AGN activity is triggered by recent infall at the cluster outskirt. In a recent work. lRuderman & Ebelinj 
(2005) study the spatial distribution of a sample of 508 X-ray point sources detected in the s oft 0.5-2.0 keV ban d in Chandra 
observations of 51 clusters (0.3<z<0.7) belonging to the MAssive Cluster Survey (MACS; lEbeling et al ■1 l200ll) . The surface 
density of the X-ray point sources computed in the cluster rest frame exhibits an 8 cr excess within 3.5 Mpc of the cluster 
centers. The authors resolve two distinct components of the excess, namely a central excess of AGN (at ^ 4 cr within 0.5 Mpc, 
as can be seen from their Fig. 2) and a broad secondary excess observed at about the virial radius, separated by a depletion 
region around 1.5 Mpc. They suggest that the central excess may be due to galaxy mergers and tidal interactions involving the 
central giant elliptical galaxies. The second excess could be caused by increased AGN activity at the cluster-field interface due 
to merger-induced accretion onto massive black holes. 

In conclusion it seems that X-ray selected AGN are broadly distributed across the cluster. This could be a consequence of 
the different processes occurring close to the center of the cluster with respect to those occurring at the cluster-field interface. 
As mentioned earlier, radio galaxies are instead very centrally peaked. These two different distributions (X -ray AGN vs radio 
galaxies) are not necessarily at variance since at the typical values of the radio luminosity (see for instance iFalcke et alJl2004l) 
the X-ray emission drops below the current X-ray detection levels. 

In order to explore and investigate further whether and how the cluster environment plays a role in the statistical AGN 
properties of galaxies, the present paper focuses on the properties of point sources detected in the inner region of high-z clusters 
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using the X-ray energy band. Our X-ray analysis is limited to the central < 1 Mpc cluster region, so as to detect point sources 
embedded in the diffuse emission of the cluster gas. Since the number of sources expected in this area is not statistically significant 
for single cluster studies, we used a combined sample of eighteen clusters observed by Chandra. The redshift range is chosen 
to cover moderate to high redshift clusters (0.25 < z < 1.01) in order to look for a possible redshift dependence of any excess 
found. All errors in this paper a re at the 1 <x confidence level, unless otherwise noted. Throughout the paper, we use a ACDM 
cosmology dSpergel et alj|2003b with H = 70 km s" 1 Mpc 1 and Q m = 1 - C2 A = 0.3. 



2. The cluster sample 

The sample used in this study consists of 18 clusters of galaxies observed by the Chandra observatory with redshift in the range 
z = 0.25-1.01. The choice of the redshift interval (z > 0.25 and up to z ~ 1) was dictated by the requirement to select clusters that 
are reasonably matched in size by the field of view of the Chandra CCDs and, at the same time, at a cosmologically significant 
distance. We retrieved only those observations from the Chandra archive with exposure times greater than 30 ks to be able to 
detect low surface-brightness clusters and the faint point sources projected against them. The choice of the exposure time was 
also dictated by the requirement of having enough count statistics to accurately measure the cluster temperatures that we will 
present in a future paper (Branchesi et al. in preparation) where the effect of the X-ray point sources on the cluster properties will 
be addressed. Table Q] lists the sample parameters and details of the Chandra observations. The columns contain the following 
information: 

- Column 1 : Cluster name 

- Column 2: Spectroscopic redshift tabulated in the literature 

- Column 3-4: Right ascension and declination (Equatorial J2000, HH MM SS.S, +DD MM SS.S) of the centroid of the 
Chandra photon distribution in the 0.5-5 keV energy band assumed as the cluster center 

- Column 5: Identification number of the observation 

- Column 6: Detector where the aimpoint lies (I, for ACIS-I or S, for ACIS-S) 

- Column 7: Observation mode (F for FAINT or VF for VFAINT) 

- Column 8: Exposure time in ks corresponding to the nominal exposure filtered to exclude time periods of high background 

- Column 9: Column density of Galactic hydrogen in units of 10 20 cm -2 , obtained from the Chandra X-ray Center (CXC) 



Propo sal Planning Tool Colden (Galactic Neutral Hydrogen Density Calculator): NRAO-compilation bv lDickev & Lockman 
dl990l) . 

- Column 10: Estimate of the luminosity limit for a cluster point source in erg s" 1 in the 0.5-2.0 keV energy band. This limit 
is computed using the flux corresponding to 90% of the sky coverage of each cluster (see Appendix[C]for details) 

- Column 1 1 : Same as Column 10 but for the 2.0-10.0 keV energy band 



With the exception of six clusters coming from the Einstein Medium Sensitivity Survey (EMSS; Gioia et al. 1990), all the clus- 
ters were originally discovered in ROSAT surveys, either the ROSAT All Sky Survey (IVog es et all 19991) or serendip i tous surveys 



from pointed observations. Four clusters come from the 160 Square Degree s ROS AT Survey (Vikhlin in et al.lll9 98: Mulli s et al. 
20031) . three from the ROSAT Deep C luster Survey (RDCS; (Rosati et al l 1 19981). on e from the Wide Angle ROS AT Pointed 



Survey (WA RPS; 



2003; Henrv et al 



Perlm anet ai1l2002l) , and one, RXJ 1716+6708 dGioia et al.lll999l) . is part of the NEP survey (Gio ia etal 



2006). Three clusters, ZWCL 0024.0+ 1652, ZWCL 1454.8+2233, Abell2125 are instead optically selected 



clusters. 



3. Source detection strategy 

Data reduction was performed using version 3.2. 1 of the CIAO software (Chandra Interactive Analysis of Observations; see web 
page http://cxc.harvard.edu/ciao/index.html). The details are given in AppendixlAl 

For each ACIS CCD chip, two separate images were extracted from the event file for source detection at the raw resolution of 
0.492 arcsec pixeF 1 . The two images are characterized by the following energies: a soft-energy image (0.5-2.0 keV) and a hard- 
energy image (2.0-7.0 keV). The cut-off below 0.5 keV is necessary due to the steep drop off of the quantum efficiency and to the 
steep rise observed in the background rate due to charged particles. The cut-off above 7.0 keV is necessary due to the decrease 
in the effective area of Chandra and to the increase in the instrumental background, which limits the detection efficiency of sky 
and source photons. 



Sources were detected using the WAVDETECT algorithm (Fre eman et al.ll2002l) . included in the CIAO software package. The 
significance threshold used for source detection was set to the inverse of the total number of pixels, e.g. ~ 10~ 6 for a 1024 x 1024 
pixels field. This is equivalent to stating that the expected number of false sources is one over the area of a single full-resolution 
ACIS chip. Wavelet scales were chosen in nine steps of ( y2)' pixel (i = 0, .., 8) starting from 0.492" , i.e. (0.5" - 8.0" ) to cover 
a wide range of source sizes, accommodating extended sources and the variation in the PSF as a function of the off-axis angular 
distance, (i.e. the distance of the source from the aimpoint). 
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Table 1. Cluster sample parameters and details of Chandra observations 
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The algorithm also uses an exposure map for each energy band to account for variations in the effective exposure across the ACIS 
field of view. To consider the photon-energy dependence of the effective exposure time, exposure maps were created at a single 
energy resolution representative of the mean energy of the photons in each band: 1.0 keV for the soft band and 4.0 keV for the 
hard band. 



4. X-ray point-source sample 

Even if WAVDETECT is a detection algorithm, it supplies a list of source parameters that is very useful for photometric analysis. 
After running the algorithm on the different energy images, these parameters were used to build a sample of point sources. 
Following Ma nners et al. I J2003I) and lJohnson eTail d2003l) . the source lists were built by accepting sources with a signal-to-noise 



ratio (S /N) greater than 3.0. This limit is a reasonably conservative one that guarantees the reliability of the sample sources. The 
source significance is defined as 



S/N = C/(l + VO.75 + B) (1) 

where C are the net source counts, and B are the b ackground counts w ithin the "source cell", a region defined by WAVDETECT, 
which is assumed to contain all the source co unts (jFreeman et al.ll2002 ). The denominator of Eq.[T]is an approximate expression 



for the error on the background counts from lGehrelsl ( 1 1986b . who gives the upper confidence level equivalent to 1 cr Gaussian 
error for small number statistics. The definition of source significance is computationally convenient for defining a flux limit 
(see Appendix in each energy band. For a number of clusters we noticed that the detection algorithm tends to consider as 
point sources some slightly extended emission regions close to the cluster center. In most cases, these sources are simply cluster 
clumps of the thermal gas rather than central X-ray point sources. A visual inspection of all the detected sources enabled us to 
tentatively eliminate these dubious point source identifications from the sample. 

A total of 119 X-ray sources were detected in the searched cluster area. Of these, 41 sources were detected only in the soft 
band, 24 only in the hard band, and 54 are common to both bands. The search radius R ext (listed in Table [2] for each cluster) is 
the radius that includes the cluster region where diffuse emission is still detectable. It is namely the radius at which the cluster 
surface-brightness profile merges into the background, and beyond which no further significant cluster emission is detectable. 
The total area covered by the clusters is ~ 0.083 deg 2 . The survey is complete down to a flux limit of 2.7 x 10~ 15 erg crrT 2 s~' 
in the soft energy band and to 0.8 x 10~ 14 erg crrT 2 s~' in the hard energy band, corresponding to 100% of the respective sky 
coverages (see AppendixO. However there are a number of individual clusters surveyed to a fainter point source flux. 
The source list is given in Table [2] A detailed description of how the X-ray source parameters are computed is given in 
AppendixlBl The columns of Table|2]contain the following information: 



- Column 1 : Cluster name 
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- Column 2: Search radius, R ext , in arcsec 

- Column 3: Source identification number 

- Column 4-5: WAVDETECT source position; Right Ascension and Declination (Equatorial J2000, HH MM SS.S, +DD MM 
SS.S) 

- Column 6-7: Net counts in the soft (0.5-2.0 keV) and hard (2.0-7.0 keV) energy bands. Asterisks indicate sources not used 
in Sect. [5] for the computation of the source counts since their flux is smaller than the flux corresponding to 20% of the sky 
coverage (see AppendixO 

- Column 8: Observed X-ray flux in units of 10~ 15 erg ctrr 2 s _I in the soft (0.5-2.0 keV) band. The flux and the associated 
error have been calculated as described in AppendixlBl 

- Column 9: Correction factor to be applied to the observed soft X-ray flux to obtain the Galactic absorption corrected X-ray 
flux. 

- Column 10: Same as Column 8 but for the hard (2.0-10.0 keV) band. 



Table 2: X-ray cluster region point sources 
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# 


RA 
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ZWCL 1454+2233 200 


1 


14 57 9.7 


+22 23 04.0 


22.10 


14.17 


1.69 + 0.40 


1.096 


4.29 ± 1.32 




2 


14 57 10.8 


+22 18 45.0 


46.94 


21.56 


3.13+0.47 




6.11 ± 1.47 




3 


14 57 12.2 


+22 21 42.4 


87.63 


27.32 


6.41 ±0.76 




8.10 ± 1.90 




4 


14 57 13.2 


+22 17 27.0 


25.65 


15.59 


1.71 ±0.34 




4.43 ± 1.23 




5 


14 57 14.8 


+22 19 33.5 


29.82 




2.07 ± 0.47 








6 


14 57 17.7 


+22 19 22.8 


35.03 


17.83 


2.58 ± 0.48 




5.46 ± 1.48 




7 


14 57 21.0 


+22 23 35.3 


88.99 


22.61 


7.44 ± 0.89 




7.15 ± 1.86 


MS 1008-1224 172 


1 


10 10 21.4 


-12 40 07.9 


5.90* 


7.84 


0.69 ± 0.29 


1.203 


4.49 ± 1.69 




2 


10 10 24.7 


-12 40 16.9 




13.74 






8.81 ±2.57 




3 


10 10 26.4 


-12 38 10.9 


7.50 




0.92 ± 0.35 








4 


10 10 29.0 


-12 40 13.5 




25.93 






15.39 ±3.47 




5 


10 10 32.3 


-12 39 34.8 


17.68 




2.21 ±0.66 








6 


10 10 35.3 


-12 40 22.0 


47.40 


16.31 


6.45 ± 0.97 




10.67 ±2.91 




7 


10 10 37.1 


-12 38 57.8 


16.82 




2.47 ± 0.66 







6 
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Table 2: X-ray cluster region point sources 



V-l LIS LCI 1 Rl 1 111. 


R ext 


# 


RA 


DEC 


^0.5-2.0 


W.0-7.0 


c 

>J0. 5-2.0 


f \ T 


c 

>J2. 0-10.0 




tt 




hh mm ss.s 








1U cgs 




1U cgs 






\->) 


V*) 


w) 
















8 


10 10 39.4 


-12 41 09.2 


9.40 




1.24 ± 0.42 










9 


10 10 42.7 


-12 39 19.1 


20.88 




2.84 ± 0.65 










10 


10 1041.8 


-12 40 02.3 


11.84 




1.59 + 0.50 






ZW CL 0024+ 1 652 


118 


1 


00 26 3 1 .0 


1 i~l 1 f\ O f\ A 

+ 17 10 30.4 


15.71 




1 o a r\ O A 

1.34 ± 0.34 


1.125 








2 


00 26 31.1 


+ 17 10 17.3 


235.11 


66.39 


20.81 ± 1.38 




43.79 ± 7.25 






3 


00 26 31.7 


+ 17 10 22.6 


16.70 




1.48+0.36 










4 


00 26 32.0 


+ 17 09 41.9 




19.60 






12.52 ± 3.17 






5 


00 26 33.0 


+ 17 07 59.9 


33.51 


6.70 


2.93 + 0.51 




4.38+ 1.80 






6 


f\f\ O/C O O o 

00 26 33.3 


■ 1 '"7 1 A O A C 

+ 17 10 34.5 




c n a 

6.74 






4.27 ± 1.74 


MS 1621+2640 


148 


1 


16 23 29.0 


+26 34 46.8 


5.81* 




1.01 ± 0.43 


1.107 








2 


16 23 33.9 


+26 35 24.9 


7.65 


7.76 


1.32 + 0.49 




6.45 ± 2.45 






3 


1 C OO A A O 

16 23 40.3 


■ o/c O C C A A 

+26 35 50.0 




1.11 






C C O i O A A 

6.58 ± 2.49 






4 


16 23 43.7 


+26 32 44.7 


282.84 


107.61 


52.80 + 3.37 




94.81 ± 13.50 






5 


16 23 45.8 


+26 33 35.2 




15.34 






13.44 ± 3.77 


RXJ 1701+6414 


108 


1 


17 01 13.1 


+64 12 50.4 


5.88* 


7.61 


0.67 ± 0.28 


1.077 


3.94 ± 1.52 






2 


17 01 21.5 


+64 15 05.2 




8.54 






4.50+ 1.65 






3 


17 01 28.3 


+64 13 32.7 




32.31 






16.54 ± 3.40 


CL 1641+4001 


89 


1 


16 41 50.0 


+40 02 49.4 


13.82 


7.79 


1.99 ± 0.54 


1.030 


4.57 ± 1.72 






2 


16 41 50.3 


+40 01 45.7 


162.94 


46.56 


22.74 ± 1.95 




26.83 ±4.83 






3 


16 41 54.2 


+40 00 32.6 


172.48 


93.56 


23.72 ± 1.98 




53.58 ±7.85 






4 


16 41 53.6 


+40 01 45.0 


8.71 




1.22 + 0.45 










5 


16 41 55.4 


+40 01 43.0 




7.59 






4.38 ± 1.70 






6 


16 41 58.4 


+40 00 48.4 


9.89 


6.89 


1.38 ± 0.44 




3.98 ± 1.58 


V 1524+0957 


148 


1 


15 24 30.6 


+09 57 30.5 


12.77 


13.76 


1.63 + 0.46 


1.087 


7.20 ± 2.09 






2 


15 24 32.3 


+09 57 45.1 


23.74 


18.80 


2.97 ± 0.62 




9.75 ± 2.48 






3 


-l c r\ A "\ ^ A 

15 24 32.4 


+09 59 07.2 




5.91 






3.39 ± 1.45 






4 


15 24 35.5 


+09 58 22.3 


7.82 


9.74 


0.99 ± 0.36 




5.24 ± 1.79 






5 


15 24 38.0 


+09 58 53.2 


21.56 


10.81 


2.55 ± 0.56 




5.45 ± 1.77 






6 


15 24 42.4 


+ 10 00 01.0 


6.89 


6.83 


0.85 + 0.33 




3.69 ± 1.48 






1 


1 C O A AO A 

15 24 43.4 


■ A A C C O C A 

+09 55 36.0 


oo 1 i ^ 

221.31 


AC TO 

95.38 


O O O C i A OO 

28.86 ± 0.22 




CA Ol , H AC 

50.83 ± 7.45 






8 


15 24 43.7 


+09 56 05.1 


7.55 




0.96 ± 0.36 






MS 0451.6-0305 


148 


1 


04 54 12.9 


-03 00 46.8 


42.92 




3.33 ± 0.58 


1.156 








2 


04 54 10.9 


-03 01 24.3 


13.26 




1.13 ±0.35 










3 


04 54 12.3 


-02 59 11.3 


6.74 


7.78 


0.54 ± 0.21 




4.47 ± 1.70 






4 


f\ A C A 1 aC C\ 

04 54 16.0 


Al AO OO O 

-03 02 32.2 




6.76 






A A 1 i 1 tl A 

4.01 ± 1.64 


V 1121+2327 


128 


1 


11 20 49.7 


+23 27 21.2 


54.77 


26.67 


4.70 ± 0.65 


1.039 


10.10 ± 2.23 






2 


11 /in O 

1 1 20 49.8 


■ O O O tl 1 A A 

+23 26 30.4 


8.83 


13.83 


a nr A O C 

0.75 ± 0.25 




5.17 ± 1.50 






3 


1 1 20 54.0 


+23 27 04.9 


160.70 


44.51 


13.98 ± 1.19 




16.89 ± 3.09 






4 


11 20 58.8 


+23 26 29.6 




8.48 






3.25 ± 1.20 






5 


11 21 04.7 


/-> ^ /-> r* -i -i A 

+23 25 11.4 




6.78 






2.64 ± 1.07 


MS 2053-0449 


118 


1 


20 56 14.3 


-04 37 16.8 


9.81 




1.36 ± 0.44 


1.149 








2 


20 56 18.7 


-04 39 14.6 


7.91 




1.03+0.37 






V 1221+4918 


143 


1 


12 21 12.6 


+49 19 19.1 


6.68 




0.52 ± 0.21 


1.043 








2 


12 21 18.1 


+49 16 35.6 




7.82 






2.47 ± 0.93 






3 


12 21 20.1 


+49 18 44.0 


40.51 


153.94 


31.03 ± 1.83 




50.08 ± 6.59 






4 


12 21 26.3 


+49 18 04.1 




12.18 






3.92 ± 1.23 






5 


12 21 30.9 


+49 17 57.3 


9.29 




0.76 ± 0.26 










6 


12 21 29.1 


+49 16 43.4 


7.82 


8.72 


0.57 + 0.21 




2.76 ± 9.92 


MS 1137+6625 


103 


1 


11 40 06.2 


+66 08 18.2 


6.99 




0.30 + 0.12 


1.035 








2 


11 40 12.8 


+66 07 33.0 


9.64 


21.33 


0.41+0.14 




4.70 ± 1.15 






3 


1 1 40 20.4 


+66 07 30.5 


39.55 


22.19 


1.67 + 0.27 




4.87 ± 1.18 






4 


11 40 31.2 


+66 08 58.2 


860.74 


263.18 


36.44+ 1.43 




59.46 ± 7.50 
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Table 2: X-ray cluster region point sources 
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V^lUSLcr Ilalllc 


D 

R ext 


a 

7T 


R A 




M).5-2.0 


^2.0-7.0 


c 

JO.5-2.0 




c 

^2.0-10.0 




tt 




hh mm ss.s 


o t If 






1U cgs 




iu cgs 


( 1 ) 


(2) 


(3i 


(4) 




(6) 


(7i 


\°j 


(9) 








5 


11 40 33.7 


+66 07 39.6 


7.71 


19.63 


0.32 ± 1.18 




4.35 ± 1.10 


RDCSJ 1317+2911 


69 


1 


13 17 18.9 


+29 11 11.1 


62.43 


24.59 


3.22 ± 0.42 


1.031 


5.54 ± 1.27 






2 


13 17 20.7 


+29 12 01.6 


6.85 




0.35 ±0.15 










3 


13 17 22.0 


+29 1 1 24.2 




21.58 






4.89 + 1.18 






4 


13 17 23.5 


+29 11 49.5 


7.70 




0.41 ± 0.15 






RDCSJ 1350+6007 


128 


1 


13 50 37.8 


+60 08 21.2 


39.73 


13.78 


4.15 ± 0.67 


1.052 


6.11 ± 1.78 






2 


13 50 39.8 


+60 05 06.8 




6.85 






2.95 ± 1.18 






3 


13 50 43.0 


+60 06 09.3 


6.92 




0.68 ± 0.26 










4 


13 50 46.1 


+60 06 58.2 




11.65 






5.03 + 1.58 






5 


13 50 50.2 


+60 08 01.5 


7.86 




0.98 + 0.35 










6 


13 50 50.4 


+60 06 20.5 


6.72 




0.67 ± 0.26 










7 


13 50 57.0 


+60 07 28.6 




9.79 






4.59 ± 1.56 






8 


13 50 57.7 


+60 08 13.7 


123.40 


21.74 


13.04+ 1.25 




9.63 ± 2.31 






9 


13 51 04.6 


+60 06 27.5 


11.89 


10.86 


1.22 ± 0.37 




4.76+ 1.54 


RXJ 1716+6708 


108 


1 


17 16 36.9 


+67 08 30.0 


8.71 


55.42 


1.05 ± 0.36 


1.111 


32.1 ± 5.48 






2 


17 16 37.6 


+67 07 31.0 


8.85 




1.02 + 0.34 










3 


17 16 42.2 


+67 06 59.8 


17.33 




1.84 + 0.45 










4 


17 1651.7 


+67 08 54.8 


55.57 


54.10 


5.63 ± 0.77 




26.71 ±4.60 






5 


17 16 53.1 


+67 07 50.2 


46.09 


16.52 


4.75 + 0.71 




8.28 ± 2.24 






6 


17 17 07.4 


+67 08 40.0 




5.85* 






2.88 ± 1.24 


MS 1054-0321 


128 


1 


10 56 51.4 


-03 38 00.7 


7.77 




0.35 ± 0.13 


1.110 








2 


10 56 52.6 


-03 38 19.8 


23.59 


7.75 


1.06 + 0.22 




2.65 ± 1.01 






3 


1 A C£ CO O 

10 56 58.8 


m o o c 1 h 
-03 38 51.2 


^Tl 1 a 

z/z.iv 


1 O A (11 

124.91 


Yl.Zy ± 0. / / 




A H AO i £. AO 

42.08 ± 6.08 






4 


10 57 04.9 


-03 38 21.2 




9.71 






3.31 ± 1.14 


WARPJ 1415+3612 


79 


1 


14 15 11.9 


+36 11 24.7 


16.59 


10.72 


1.11 + 0.28 


1.033 


3.03 ± 0.99 






2 


14 15 12.4 


+36 13 03.9 


8.70 




0.60 + 0.21 










3 


14 15 13.5 


+36 12 10.4 




22.51 






6.39 ± 1.52 






4 


14 15 16.1 


+36 11 51.8 


40.76 


19.70 


2.77 ± 0.45 




5.62 ± 1.40 



5. X-ray point-source number counts 

We used the source counts, or log N- log S relationship, as a statistical tool for investigating whether there is any point source 
count excess in the regions of diffuse cluster emission with respect to the fields void of visible clusters. Although integral source 
counts are usually given in the literature, we present here both the differential and integral distributions. The differential number 
counts are statistically more correct since source-count errors are independent, di fferently from the integral number counts. To 



compute the log TY-log S relationship, we followed the method described in iGioia et al.l {[990). The source contributions 
are computed by properly weighting each source with its visibility area, i.e. the area of sky where the source intensity equals 
or exceeds the sensitivity limit (see Appendix ICl for details). The integral log (N > S) diagram is then built by summing, in 
decreasing flux order, the contribution from each source. Since each cluster sky coverage rapidly decreases near the limiting flux 
for both the soft and hard bands, in order to prevent incompleteness effects, we considered only those sources with a flux larger 
than the flux corresponding to 20% of the total sky area covered by each cluster. It is worth noting that, throughout the following 
sections, the soft X-ray fluxes are unabsorbed fluxes, i.e. corrected for the Galactic hydrogen column density absorption along 
the line of sight at each field position. 



5.1. X-ray source counts in the cluster regions 

Following lCrawford et al. dl970h and lMurdoch et alj dl973l) , a Maximum Likelihood Method (MLM) that operates on the differ- 
ential counts, was adopted to determine the slope of the number-flux distribution of the X-ray sources. The first assumption is 
that the differential log N- log S distribution may be described by a single power-law model of the form: 
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The MLM uses the 'unbinned' data that is very useful when one deals with a low number of sources. 

The MLM applied to the sources found in the regions covered by our clusters (from now on cluster region sources) provides the 
following results: 

dN I S \-178±0.08 

dS = 354 ± 37 bxlO^ergcm-^ ) S ™ ^ (3) 

for the 0.5-2.0 keV band and 

dN I S \-2.0±0.1 

— = 255 + 29 —-. -. — r sources deg" 2 (4) 

dS \1 x 10- 14 erg cm- 2 s- 1 / 5 

for the 2.0-10.0 keV band. 

The values for b were calculated at the flux levels of So — 2 x 10~ 15 erg crrT 2 s~' in the soft band and So = 1 X 10~ 14 erg cm~ 2 
s~' in the hard band. For the exposures considered here, these values correspond to the centers of the sampled flux intervals. 

The differential logarithmic source counts, in bins of A log S = 0.2, are indicated in Fig. [T]by filled triangles. The soft (hard) 
log N- log S for the X-ray point sources detected in the areas covered by the clusters is given in the top (bottom) panel. The 
counts and the errors are given by: 



n j n 

dN=Y— deg 2 and cr dN = « Y(l /Q 2 ) deg 

i=i i2 ' > i=t 



(5) 



where n is the number of sources within each bin and O, the visibility area of the ith source. The thick solid lines represent the 
best fit obtained with the MLM described above. The dashed lines indicate the 1 cr uncertainties on the MLM fit parameters. 

To compare our results with those in the literature, we also calculated the cumulative source number counts N (> S) and the 
corresponding errors as follows: 



N(> S) = V — deg" 2 and o- N(>s) = . V (1/Q 2 ) deg" 2 

i=i lli \ t=i 



(6) 



where N m is the total number of detected sources with a flux > S and Q, the sky coverage corresponding to the flux of the ith 
source. The assumption made of a single power law converts Eq.|2]into 



N(>S) = k \S-o) (?) 

for the integral distribution, where or = al - 1, The integral normalization k is calculated by integrating Eq.|2]between the flux 
limit and infinity: 

The errors on k were calculated using the propagation of error for al and b. The results for the integral source counts in the 
cluster fields are listed in Table [3] 

Figure [2] illustrates the cumulative source counts, indicated by solid triangles, which were plotted using a step of A log S = 
0. 15. The soft integral log N— log S is shown in the top panel and the hard integral log N— log S in the bottom panel. Errors 
correspond to 1 cr confidence level and are obtained using Eq. [6] Note that the errors on the integral log N— log S are not 
independent. 



5.2. X-ray source counts in the reference fields 

The slope and/or normalization of the source counts at a given flux limit can provide information on the point-source population 
in the cluster regions compared to the sources in fields without visible clusters (from now on field sources). To check for the 
presence of such an effect, the results obtained in the area of the 18 clusters were compared with those obtained in four reference 
fields void of clusters with exposure times similar to the exposures of the cluster fields. 

The reference fields used were the Chandra Deep Field South (CDFS), the Hubble Deep Field North (HDFN), the Groth/Westphal 
strip area, and the Bootes field (see TableO. In order to sample similar flux levels, the analysis was limited to observations with 
exposure on the order of lOOks. The analysis was performed on the four ACIS-I CCD (16.8'xl6.8')- For the validity of our 
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analysis and to avoid systematic errors, we ran exactly the same procedure on the reference fields as was used to derive the 
log N— log S of the cluster region sources. The same source-detection routine, flux limit and flux estimate method have been 
adopted (see Sects. [3] and [4] and Appendix for details). The sky coverage of each reference field was constructed following the 
method described in Appendix|C] The only difference is the binning of the background and exposure maps, which were chosen 
to be 64 x 64 pixels to reduce the computer time. The value of the flux limit of each bin is a representative mean of all the original 
pixels since the absence of a galaxy cluster in the field results in a much more uniform background. The sky area of each field 
measures ~ 0.078 deg 2 . Following the same approach as the cluster region sources we considered only field sources with a flux 
greater than the flux corresponding to 20% of the sky coverage of each field. 



The best-fit results on the cumulative log N— log S in each reference field are summarized in Table [3] Note that the literature 
reports source counts that are usually estimated using the observed fluxes, while we corrected for the Galactic hydrogen column 
density to produce unabsorbed fluxes. The counts agree with recent studies of these fields dGiacconi et al. 2001 ; Tozzi et al 



Table 3. Integral source counts in clusters and in reference fields 



Name 


Obs.Id 


r 


NT* 


No.5-2.0 


N2.0-7.0 


co.5-2.0 


o^.o-io.o 


ko.5-: 


.0 


ki.o-io.o 






ks 


cirr 2 










deg- 


2 


deg" 2 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 




(10) 


Clusters 








92 


77 


0.78 ± 0.08 


1.00 ± 0.10 


454 ± 


66 


255 ± 39 


Bootes 


3130 


107 


1.15 


90 


84 


0.76 ± 0.06 


1.11 ± 0.10 


410 ± 


54 


217 ± 31 


HDFN 


3389 


109 


1.50 


96 


63 


0.78 ± 0.08 


0.90 ± 0.10 


426 ± 


62 


240 ± 41 


CDFS 


582 


128 


0.80 


109 


81 


0.80 ± 0.06 


1.03 ± 0.10 


396 ± 


48 


210 ± 31 


Groth 


4357 


85 


1.30 


83 


52 


0.80 ± 0.08 


0.94 ± 0.12 


432 ± 


63 


217 ±41 


Unified 








378 


280 


0.79 ± 0.04 


1.01 ± 0.04 


416 ± 


30 


219 ± 16 



Column 1 : Reference field name 

Column 2: Identification number of the observation 

Column 3: Total net exposure time after 'flare' cleaning 

Column 4: Hydrogen Galactic column density value, N#, in units of 10 20 cirr 2 . 

Column 5-6: Number of sources detected in the soft and hard bands with flux larger than the flux corresponding to 20% of the sky coverage 
Column 7-8: The integral log N— log S power law slope and 1 cr confidence level error for the soft and hard bands 

Column 9-10: The integral log N- log S normalization and 1 cr confidence level error calculated in the soft band at a flux of 2 x 10~ 15 



erg cm s and in the hard band at a flux of 1 x 10 
Galactic hydrogen column density. 



erg cm 2 s . Note that the soft band normalization takes into account the different 



200ltlRosati etalj2002 ; Bran dt et al 1 12001 ; Wan g et al . 2004). The slope values match the published ones within the errors. The 
normalizations are also consistent within the errors, considering that we used different, more recent calibrations and a different 
spectral slope (T = 1.7) to derive the conversion factors (see Appendix IBl. Most of the authors quoted above use a flatter slope 
with r = 1.4. The average difference between the conversion factors obtained considering an absorbed power law spectrum with 
r= 1.7 andT = 1.4 are 



K soft (lA) 



-tf„ /t (1.7)- 



- 3%K soft {l.l) 
\2%K hard (\.l). 



Finally the four reference fields were combined in a unique big field of ~ 0.3 1 deg 2 . A corresponding sky coverage was also con- 
structed. The 378 plus 280 sources, detected in the soft and hard band respectively, have been used to derive both the differential 
and integral log N— log S for the unified reference field. The data were best-fitted using the MLM. The best-fit integral source 
count parameters are listed in Table [3] 

Hereafter the unified reference field is used for comparison with the cluster fields, but the four single reference fields give a 
measure of the parameter dispersion. The differential source counts in bins of A log S = 0.2 for the unified reference field are 
indicated by open squares in Fig. Q] The errors correspond to 1 cr confidence level obtained as in Eq. [5] The thin solid lines 
represent the power law corresponding to the best-fit parameters obtained with the MLM. 

The cumulative source counts of the unified reference field are shown in Fig. [2] (open squares). The dashed lines represent 1 cr 
uncertainties on the source counts. 
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6. Are X-ray point sources associated with the clusters? 

6.1. Cluster counts vs field counts 

Figure [3] illustrates the best-fit normalization and slope of the integral log N- log 5 s, for both the cluster and the reference field 
sources. The parameters of the fits, given in Table [3] were obtained by fitting one single power law over the entire flux range. 
All values are consistent within 1 cr confidence level in both bands, showing no significant difference between the counts in 
the cluster areas and in the control fields. However, an inspection of the top panel of Fig. [2] shows that the soft integral cluster 
counts exhibit a small excess at the bright end around ~ 2 x 10" 14 erg crrT 2 s -1 . A similar excess is present in the hard integral 
number counts (Fig.|2j bottom) around ~ 3.0 x 1CT 14 erg cm -2 s _I . The significance of these apparent excesses was evaluated by 
computing the number of field point sources expected in the cluster regions. Since the fie ld source counts in both bands show a 



steepe r slope at the brigh t end compared to the faint end, as also found by several authors (IHasinger et alJ I 1998: Mus hotzky et al 



2000; iTozzi et al]l2001t iRosati et al] 120021) . instead of using the single power-law best fit, we considered it more appropriate 



computing the expected number of field sources by directly normalizing the number of the observed point sources in the control 
fields to the area covered by the clusters. We would expect 5+1 field sources above a flux of 1.25 x 10~ 14 erg cm -2 s _1 for the 
soft band and the same number above 2.5 x 10~ 14 erg crrT 2 s _1 for the hard band. We instead find 1 1 sources in the cluster region 
in each band. The significance of the difference is at the 1.7 cr level in both cases. Note that 9 sources out of the 1 1 are common 
to both bands. If this small excess is real, we are facing a cluster population of X-ray sources with X-ray luminosities ranging 
from « 3 x 10 42 to 10 44 erg s _1 in the soft band and from * 10 43 to 10 44 erg s~ l in the hard band. 

Next we examined whether there is any dependence of the number counts on cluster redshift. For this purpose the point source 
sample was divided into two subsamples: the "high-z" subsample corresponding to the sources detected in the regions covered 
by clusters with redshift z > 0.7 and the "low-z" subsample corresponding to the sources detected in the regions covered by 
clusters with redshift z < 0.7. We performed an analysis on these two subsamples similar to the one performed on the whole 
sample. No significant difference between the two subsamples is evident in the soft band. In the hard band there is instead a slight 
indication of a possible excess at bright fluxes in the "high-z" subsample with respect to the "low-z" one. The statistics are too 
poor to draw any conclusion. 

Finally it is interesting to compare our estimate of the "cluster X-ray population" with the results by lMartini et al. I d2006h for a 



sample of nearby clusters. Under the assumption that our clusters are similar to the Martini clusters, we used their spectroscopi- 
cally confirmed cluster X-ray sources to derive a "Cluster X-ray Luminosity Func tion" (Cl-XLF) and predict how many cluster 
X-ray sources we should expect in our sample. In their five clusters with z > 0. 15, Mar tini et al] ( 20061) list seven objects within 
1 Mpc from the cluster center with a soft L x > 10 42 erg while no object is present in the three closer clusters with redshifts 
of 0.07, 0.06, and 0.059. If we build the Cl-XLF from the five z>0.15 clusters by properly accounting for our luminosity limits, 
we expect five cluster X-ray sources, four of which are within 0.5 Mpc and one between 0.5 and 1 Mpc from the cluster centetQ. 
These expected numbers agree with the numbers of sources in clusters previously estimated from our own data. Had we used 
the data from all the Martini eight clusters, the Cl-XLF would have had a lower amplitude, and the predicted number of cluster 
X-ray sources in our sample would ha ve dropped to < 4. A similar analysis cannot be performed in the hard band due to the 
presence of only two X-ray sources in Marti ni et al. (2006) work. 



6.2. Radial distribution of the point sources 

To further investigate the reality of the excess at high flux, we examined the radial distribution of the brightest sources in both 
bands as a function of the projected linear distance from the cluster center. Given the small number of sources, only two 0.5 
Mpc wide annuli were used. From the control field source counts, normalizing to the areas actually examined, we computed the 
expected number of field sources in the two cluster annuli. The results are listed in Table [4] (Columns (2) and (5)). In each band 5 
field sources are expected in the area covered by the clusters (see Sect. 16. II ). Two field sources are expected in the inner annulus 
against 7 observed, and 3 field sources are expected in the external annulus against 4 observed. There is a 2 cr excess within 0.5 
Mpc from the cluster center, thereby confirming the small excess of sources found in Sect. 16. II and providing additional support 
for the possible existence of X-ray sources belonging to the cluster (on average one every three clusters). 

Next we constructed the radial profile of the surface density, If 1 , within 1 Mpc (Columns (3) and (6) of Table [H, assuming the 
brightest sources to be at the cluster redshift and taking the areas effectively surveyed into account. The field source number 
density per Mpc 2 , Y) lekl (Columns (4) and (7)), was then derived assuming the same redshift distribution as in our cluster sample. 

The radial profile of the soft source surface density, If 1 , was compared with the distribution of lRuderman & Ebelind d2005l) in a 
subsample of the MACS clusters. The MACS survey was built to find the "most" massive clusters ever (see Ebeling et al.ll20o"lb . 



1 Since the search radius R ext does not extend up to 1.0 Mpc for all the clusters, the number of sources expected in each annulus was 
normalized by the area actually examined. 
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Table 4. Radial source distribution for the brightest sources 



R 

(1) 


No.5-2.0 

(2) 


^0.5-2.0 

(3) 


yheld 
"0.5-2,0 

(4) 


N2.0-10.0 
(5) 


"2.0-10.0 
(6) 


y field 
"2.0-10.0 

(7) 


0-0.5 


7(2) 


0.50 ±0.19 


0.15 ±0.03 


7(2) 


0.50 ± 0.19 


0.13 ± 0.03 


0.5-1.0 


4(3) 


0.18 + 0.09 


0.15 ±0.03 


4(3) 


0.18 ±0.09 


0.13 ±0.03 



- Column 1: Size of the annulus in Mpc 

- Column 2: Number of sources detected in the annulus with an unabsorbed flux brighter than 1.25 x 10~ 14 erg cnr 2 s -1 in the soft energy 
band. In parentheses the expected number of field sources over the cluster regions is given 

- Column 3: Soft X-ray surface density for the cluster region sources in Mpc~ 2 

- Column 4: Estimated soft X-ray surface density for the control field sources in Mpc~ 2 

- Column 5: Number of sources detected in the circular annulus with a flux brighter than 2.5 x 10~ 14 erg cirr 2 s -1 in the hard energy band. 
In parentheses the expected number of field sources over the cluster regions is given 

- Column 6: Hard X-ray surface density for the cluster regions sources in Mpc~ 2 

- Column 7: Estimated hard X-ray surface density for the control field sources in Mpc~ 2 



The median MACS X-ray luminosity!^ is 2.5 times our median luminosity with only 3 (out of 18) of our clusters more luminous 
than the MACS median luminosity. The flux limit we adopted in the soft band is the same as the one used by them (1 .25 x 10" 14 
erg cirr 2 s _1 ). Within 0.5 Mpc, where the excess is found, the two densities are in very good agreement (0.50 source s per Mpc 2 



from Table[4]vs 0.53 sources per Mpc 2 from their Fig. 2). Within 0.5 Mpc the excess of Ruderman & Ebelind (2005) is at w 4 cr 
with respect to their point-source density at the cluster field edges (4-7 Mpc) that they assume to be the background point-source 
densityB Between 0.5 and 1 Mpc, the two profiles are consistent within the errors. We stress t hat the field surface density obtained 
from our control fields agree very well with the one derived by Rud erman & Ebelin g (2005) from their cluster field edges. 



In order to increase the statistics, we extended our radial analysis to fainter sources. We used 31 sources with a soft unabsorbed 
flux brighter than 2.7 x 10~ 15 erg cm -2 s _1 and 25 sources with a hard flux brighter than 0.8 x 10~ 14 erg cm -2 s _1 since all clusters 
are complete down to these flux limits. A small source excess is found within 0.5 Mpc from the clusters center, which however 
does not improve significantly the previous result, indicating that cluster X-ray sources are confined to the highest fluxes. 



6.3. Optical counterparts 

We checked the literature for any optical counterpart associated with the 13 X-ray sources with a soft flux brighter than 1.25 x 
10~ 14 erg cirr 2 s _1 and/or a hard flux brighter than 2.5 x 10~ 14 erg cm -2 s . The results of this search are reported in Table[5] No 
optical data are available for four sources. Two sources are identified with cluster members with spectroscopic redshifts. Three 
are identified with optically faint objects with no redshift measurement. Finally three sources are background objects having a 
significantly higher redshift than those of the corresponding clusters. A fourth one may also be a background object on the basis 
of a photometric redshift. Clearly more optical work is needed in order to know the cluster membership of these objects. 



7. Comparison with radio sources in X-ray selected clusters 

As mentioned in the introduction, radio galaxies in both nearby and distant clusters have a centrall y peaked distrib ution. 
Branchesi et al. (2006) analyzed a sample of VLA radio sources detected in 18 X-ray selected clusters (IGioia et al J 1200 lb ex- 
tracted from the NEP survey with redshift and luminosity distributions similar to the pres ent sample. They foun d a pronounced 
peak of radio sources within 0.2 Abell radii, corresponding (in the cosmology adopted in lBranchesi et alJ l2006) roughly to the 
0.5 Mpc size of the first bin in Table |4]here. The radio source surface density associated to this peak is * 10 Mpc -2 , i.e. ~ 20 
times higher than the X-ray point source surface dens ity found here. How ever, the larger radio excess could be a selection effect 
due to the better sensitivity of the radio observations. iFalcke et alj (120041) discuss a correlation between the 3-9 keV luminosity 



and the 5 GHz core radio luminosity (vL„) in radio-loud AGN. To estimate the X-ray luminosity expected for the Bra nchesi et al 



(2006) radio galaxies, we first conv erted the total rad i o lum inosity to the core radio luminosity using the correlation between the 
core and radio powers published bv lGiovannini et al.l (11988). From the resulting core radio luminosity and the Falcke correlation, 
we find that the estimated X-ray luminosity of the NEP radio sources is about two orders of magnitude lower than the X-ray 



2 The MACS median luminosity was estimated using Fig. 10 in Ebeling et al. (2001) after converting their 0.1-2.4 keV i?0SAr luminosity 
into bolometric luminosity with their sa me assumptions on t he temperature, and taking into account the different cosmologies assumed. 

3 The significance of 8.0 cr claimed by Rud erman & Eb eling ( 2005) is actually for the area within 3.5 Mpc from the cluster centers. 
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Table 5. Optical parameters for the 13 X-ray brightest sources 



Cluster name 


# 


RA 


DEC 


ARA(x-o) 


ADec(x-o) 


z 


Notes 






hh mm ss.s 


O 1 II 


sec 








(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


Abell2125 


5 


15 40 56.4 


+66 16 28.7 


+0.0 


-0.3 


1.012 


a 


ZWCL 0024+ 1652 


2 


00 26 31.1 


+ 17 10 17.3 


+0.0 


+0.7 


0.400* 


b 


MS 1621+2640 


4 


16 23 43.7 


+26 32 44.7 


-0.1 


-0.3 





c 


CL 1641+4001 


2 


16 41 50.3 


+40 01 45.7 


— 


— 





- 




3 


16 41 54.2 


+40 00 32.6 


+0.0 


+0.6 


1.003 


d 


V 1524+0957 


7 


15 24 43.4 


+09 55 36.0 










VI 121+2327 


3 


1 1 20 54.0 


+23 27 04.9 








~ 


V 1221+4918 


3 


12 21 20.1 


+49 18 44.0 


— 


— 





- 


MS 1137+6625 


4 


11 4031.2 


+66 08 58.2 


+0.0 


+0.0 


1.269 


c 


RDCSJ 1350+6007 


8 


13 50 57.7 


+60 08 13.7 


+0.1 


-2.2 




c 


RXJ 1716+6708 


1 


17 16 36.9 


+67 08 30.0 


+0.0 


+ 1.0 


0.795* 


f 




4 


17 16 51.7 


+67 08 54.8 


-0.0 


-1.1 




g 


MS 1054-0321 


3 


10 56 58.8 


-03 38 51.2 


-0.1 


-0.2 


1.200'''' 


h 



- Column 1 : Cluster name 

- Column 2: Source identification number as in Table[2] 

- Column 3-4: X-ray source position; Right Ascension and Declination (Equatorial J2000, HH MM SS.S, +DD MM SS) as in Table|2] 

- Column 5-6: Offset between the position of the X-ray source and that of the optical counterpart (A RA in seconds of time and A DEC in 
arcsec) 

- Column 7: Spectroscopic or ph otometric (indicated by ph ) redshift. Asterisk in dicates cluster memb ership 

- Column 8: Literature source: a)lMilleret al.N2004h . b) ICzoske etal] ( feOOlh . cl lEckart et at] J2005I). d) from the Sloan Digital Sky Survey, 
e) Silverman et al. (2005), fl lGioia et al.l jl999h , g) from a Keck image taken by I. Gioia and h) dForster Schreiber et alj2006l) . 



luminosity reached by the present Chandra observations. Thus the X-ray instruments do not seem to have the sensitivity required 
to detect such faint X-ray counterparts of the radio-loud AGN. Similar arguments apply if the weaker NEP radio sources are 
sta r-forming galaxies d ue to the tight linear relations between the X-ray, radio, and far infrared luminosities found, among others, 
bv lRanametai](l2003h . 

On the other hand, the dozen X-ray sources at the bright end of the log N— log S have a high enough X-ray luminosity to 
statistically expect that radio emission f rom some of them w ould be detectable in surveys like the FIRST (Be cker et al.ll 19951) 
or the NVSS (NRAO VLA Sky Survey. ICondon et alJI 1998b . Indeed two such X-ray sources have associated radio emission. 
Source # 3 in CL 1641+4001 has a FIRST counterpart at the position 16 /, 41 m 54.24", +40°00'32.0" with a flux S 1AGHz = 5.06 
mJy. Source # 1 in RXJ 1716.4+6708 has a NVSS counterpart at the position 17 /, 16"'37.14 s , 67°08'28.8" with a flux S IA ch z = 
332.0 mJy (4C +67.26). Given the low statistics no conclusion is drawn here. 

For completeness the coordinates of all sources listed in Table [2] were cross-correlated with those of the FIRST (or NVSS 
catalogs when no FIRST data are available). Only one additional coincidence was found for source # 2 in ZWCL 1454.8+2233. 
The corresponding FIRST radio source at 14 /l 57 m 10.82 s , +22°18'44.9" has a flux of S 1A ch z = 4.87 mJy. 



8. Summary and conclusions 

In this paper we have presented an analysis of the X-ray point sources detected in the inner 1 Mpc region of 18 high-z (0.25 < z 
< 1 .01) galaxy clusters retrieved from the Chandra archive. Unlike most of the previous studies that analyzed the whole Chandra 
field around the clusters, we considered only the point sources embedded in the cluster emission, i.e. belonging to the clusters or 
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in projection. We find a small excess for the cluster sources at fluxes brighter than 1.25 x 1CT 14 erg cirT 2 s in the soft energy 
band and brighter than 2.5 x 10~ 14 erg crrT 2 s in the hard energy band. The significance of the excess is at the 1.7 cr level in each 
band (Sect. loTTl ). To further investigate the reality of the X-ray point source excess, we examined the source radial distribution as 
a function of the projected linear distance from the cluster center. A 2 cr excess was found within 0.5 Mpc providing additional 
support for the existence of X-ray sources belonging to the cluster. The excess is given by m 6 out of 1 1 sources in the luminosity 
range 3xl0 42 - 10 44 erg s 1 in the soft and 10 43 - 10 44 erg s 1 in the hard energy band. These results agree with those of previous 
studies that have detected excesses in cluster fields. 



5 x 10 45 erg s" 1 



with a median value of 



(only three clusters have a luminosity higher than 2 x 10 45 erg s '). Ruderma n & Ebelingl (2005) examined a 



The galaxy clusters presented here are in the bolometric luminosity range 10 44 
0.9xl0 45 erg s" 1 

sample of 51 clusters in a similar r edshift range (0-3< z < 0.7), but using more massive clusters than ours. They found a soft 
excess within 1 Mpc similar to ours. lJeltema et al. I (l2006l) analyzed six groups of galaxies in the redshift range 0.2 < z < 0.6 and 
obtained a 2 cr overdensity result for these lower-mass systems. There fore our results provide furth er evidence of the presence 
of a population of AGN in system s with very different masses. While iRuderman & Ebelind d2005l) did not analyze sources in 
the hard band. IJeltema et al.l ([2006) found no significant excess in the number of hard sources. Thus our study provides for the 
first time some evidence of source overdensity within 1 Mpc in both energy bands. In addition, there is a slight indication that 
the hard band excess increases with redshift, even though the small statistics do not allow us to draw any conclusion. We note, 
however, that this effect goes in the same direction as the appare nt correla tion between the am plitude of the overdensity and the 
cluster redshift found for the first time by Cappelluti et al. (2005) (see also Martini et alJl2006T) . 

Radio galaxies in both nearby and distan t clusters have a centra lly peaked distribution with surface density within 0.5 Mpc on 
the order of 10 Mpc -2 . A recent work by IBranchesi et aL ( 20061) on a sample of distant X-ray selected clusters (0.3 < z < 0.8) 
finds an excess in the radio surface density within 0.5 Mpc from the cluster center that is much higher (by a factor ~ 20) than 
the present X-ray source overdensity. Even if an excess of sources is present at both wavelengths, the much smaller amplitude of 
the X-ray overdensity could be explained by the better sensitivity of the radio observations. Higher sensitivity (and resolution) 
X-ray telescopes could provide more information on the nature of the population of AGN and/or star-forming galaxies in high-z 
clusters, which is responsible for the more pronounced excess in the radio domain. 
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Fig. 1. Differential log N- log Ss. The differential log N- log 5 of the X-ray point sources detected in areas covered by the 
clusters are indicated by solid triangles (soft band in the top panel and hard band in the bottom panel). Open squares indicate 
the differential log N- log S of the X-ray point sources detected in the field obtained by combining the four deep (r ~ lOOLs) 
fields without clusters (Bootes, CDFS, HDFN, and Groth Strip). The uncertainty on each point corresponds to a 1 cr confidence 
level. No cluster region sources are observed in the hard faintest flux bin. The arrow indicates an upper limit corresponding to 
three sources. According to Poisson statistics, we have a 5% probability of observing zero sources when three are expected. The 
thick (thin) solid lines represent the best-fit power law of the log N- log S for the cluster region (field) sources obtained with 
the MLM, while the dashed lines correspond to the 1 cr uncertainties on the fit for the cluster region sources. The uncertainties 
on the field source fit are smaller and have been omitted for clarity. 
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Fig. 2. Integral log A^-log 5 s. The integral log iV-log S for the sources detected in the regions occupied by clusters are 
indicated by filled triangles in both the soft (top panel) and hard (bottom panel) energy band. Open squares indicate the integral 
log N- log S for the sources detected in the unified reference field obtained by combining the four ~ 100L? long fields without 
visible clusters (Bootes, CDFS, HDFN, and Groth Strip). The cumulative source counts are plotted using a step A log S = 0.15. 
Uncertainties on each log N- log S point correspond to a 1 cr confidence level obtained as described in the text. Dashed lines 
represent 1 cr confidence level for the field sources counts. 
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Fig. 3. Integral log N- log S normalization vs slope for both cluster region and field sources obtained using the MLM. Solid 
triangles indicate the best-fit results for the cluster region sources (soft band to the left and hard band to the right), while the thick 
open squares indicate the unified reference field results. The thin open squares indicate the best-fit results for each of the four 
reference fields. Uncertainties correspond to a 1 cr confidence level on the fit parameters. 
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Appendix A: Data reduction 



We retrieved the level=l event files from the archive and applied t he standard processing. The CIAO tool acis process. events 
was used to apply the correction for charge transfer inefficiency (CTI Townsle y et alj|2000t Grant et al]|2005 ) and recomputation 
of event grades. To compute calibrated photon energies, acis process. events was also used to update the Advanced CCD Imaging 
Spectrometer (ACIS) gain maps with the latest version provided within CALDB (ver. 3.0.3) and to correct for its time dependence 
(T_GAIN correction). The CTI and T_GAIN corrections were applied for those chips and period observations for which they were 
available. Most of the observations were telemetered in VFAINT mode, which provides a better rejection of the particle-induced 
background. For these observations we ran acis .process. events with the option 'check_vf_pha=yes', to flag probable background 
events. After running acis process. events, the events were filtered to include only the standard event grades 0, 2, 3, 4, and 6 and 
status bits set to 0. In this way we removed photons detected in bad CCD columns and bad pixels, 'problem' events such as 
cosmic ray afterglows, and also those with bad ASCA grades (1,5, and 7). 

The final step was to examine the background light curves during each observation to detect and remove the flaring episodes. 
After excluding the point sources and cluster emission from the event file, the script Icclean was used to extract and bin the 
light curve and to calculate the average count rate. The flare detection was performed following the recommendations on the 
energy band and the bin size given in the Markevitch COOKBOOK0 We excluded those time periods when the count rate was 
not within 20% of the quiescent rate. The final count rate in the source-free regions of each observation was compared with the 
background value tabulated by Markevitch. We obtained values consistent within less than 10%. Compared to the ACIS front-side 
illuminated (FI) chips the back-side illuminated (BI) chips (SI and S3) have a higher sensitivity at low energies. However, this 
low-energy s ensitivity makes the chip more sensitive to particle events, which resu lts in more frequent background flares than in 
the FI chips (IPlucinsky & Viranill2000l Mar kevitch 1 1200 1 : Mar kevitch et alJl2003h . When the source lay in the S3 chip, the SI 
chip was accurately examined to exclude completely the flare-affected period. 

The released calibration underestimates the effective area of the Chandr a mirror by 10% j ust above the Ir M edge, probably 
because the mirror surface is contamin ated by a thin hydroca rbon layer (M arshall et al]l2003l) . To correct the effective area, a 
"positive absorption edge" described by lVikhlinin et alJ (2005) was used in the analysis of the spectra of point sources with the 
X-Ray Spectral Fitting Package (XSPEC). 



Appendix B: Computation of Source Parameters 

The net counts, C, calculated as the sum of all counts in the "source cell" subtracted by the sum of the background counts, B, 
were used to estimate the source flux. Since the majority of the detected sources have poor statistics (less than 50 total net counts 
in the soft or the hard band), an estimate of the source flux through a fit to the data with an absorbed power law is not always 
possible. Two separate conversion factors, one for each band, were then calculated to derive the flux {S ) from the observed count 
rate. 

The net count rate for each source was computed by dividing the net counts by the effective exposure time at each source position. 
For each source, the effective exposure time is given by the observation exposure time (corrected for the flares' time periods) 
multiplied by the ratio of the exposure map, averaged within the extraction region for each source to the value of the exposure 
map at the aimpoint. The vignetting correction (V) to be applied to the net counts is given by the ratio of the value of the exposure 
map at the aimpoint (expmap jiimpoint) to the value of the exposure map at the source position (expmap. source). The correction 
is done separately for each band using the exposure maps computed at energies of 1 .0 ke V (soft band) and of 4.0 ke V (hard band). 
In this way the source count rate corresponds to the count rate that the source would have if it were observed at the aimpoint. 
The soft and hard conversion factors (K) from counts (cts s _1 ) to X-ray fluxes ( erg crrr 2 s~') were derived at the aimpoint using 
the response matrices of the detector a t this position. Such conversion factors w ere computed assuming an absorbed power-law 
spectrum with a photon index F =1.7 ( Mus hotzkvll 1 984 : Wilkes & Elvis|[l987l) and assuming for the hydrogen column density, 



N//, the Galactic value along the line of sight at the source position (see TableQ] Column 9, and Table[3] Column 4). 
The soft and hard X-ray source fluxes were calculated as 

s = C x ex pm gp*impomt xK= C xyxK 
t expmap. source r 

where C are the net counts, t the flare-corrected exposure time, V the vignetting correction factor, and K is the conversion factor 
from counts s to X-ray fluxes in erg crrr 2 s appropriate for each energy band. The conversion factors allowed us to convert 
the 0.5-2.0 keV band count rate to the observed X-ray fluxes in the same band and the 2.0-7.0 keV band count rate to the 
observed X-ray flux in the 2.0-10.0 keV band. 



4 See http://cxc.harvard.edu/contrib/maxim/acisbg/data/README and 
http://cxc.harvard.edu/contrib/maxim/acisbg/COOKBOOK 
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Flux uncertainties were estimated taking into account the error on the net counts (cr NC ) computed as the square root of the total 
observed counts in the "source cell", i.e. VC + B, and the error on the conversion factor, cr K , due to the assumed different power- 
law models. The uncertainties on the conversion factor reflect the range of possible values for the effective photon index: F = 
1.4-2.0. The error on the flux then becomes: 




(B.2) 



The soft observed (i.e. absorbed) X-ray flux used to derived the source counts were multiplied by the factor, (Table[2] Column 

3) to obtain the Galactic unabsorbed X-ray flux. The correction factors were derived taking into account the Galactic hydrogen 
column density, Nh, along the line of sight of each cluster (TableQ] Column 9) and of each void of cluster field (Table[3] Column 

4) . 



Appendix C: Sky Coverage 

In Sect. [5] we derive the source counts (or log N— log S) for the point sources detected in the cluster regions. An important 
ingredient for the log N— log S is the determination of the sky coverage. The "sky coverage" is the area of sky sensitive down 
to a given flux limit as a function of the flux density. To estimate the sky coverage, Q, we constructed a flux limit map, indicating 
the flux of the faintest source that would have been included in our source list at each position of the ACIS chip within the region 
of the cluster. The flux limit map has been constructed to account for the following effects: 

- instrumental effects, such as vignetting, or increase of the point spread function size with off-axis angle. The sensitivity of 
the ACIS detectors varies significantly across the field of view; 

- background effects. The background considered also includes the extended emission from the cluster's hot gas, since we 
detected sources within the region occupied by the clusters. 

Following Eq. IB. li the sensitivity limit (S n, n in erg ctrr 2 s _1 ) at each detector position was defined as: 

S, im = — xV xK (C.l) 

T 

where Q, m are the net counts derived as in Eq. |C.2| 

For the chosen signal-to-noise ratio (S/N = 3; Section|4| one obtains from Eq.Q] 



Q im = 3 X (1 + ^0.75 + B avg x A) (C.2) 

where B avg is the average background counts per pixel, A the area expressed as the number of pixels covered by the source cell, 
and B avg xA is the local background counts within the source region. The estimate of this area A takes into account the degradation 
of the Chandra PSF, which increases with the off-axis angle. To estimate the local background (inclusive of the cluster extended 
emission) and the vignetting factor at each position of the chip, we built both a background and an exposure map of the same 
size and resolution for both the soft and the hard energy bands. The map of the effective background plus the extended cluster 
emission were obtained as follows: 

1 . All identified point sources were subtracted from the soft and hard images. The resulting "holes" were filled with pixel values 
sampled from the Poisson distribution whose mean and standard deviation equalled that of the surrounding background pixels 
(using the CI AO tool dmfilth); 

2. These source-free images were binned by 32 x 32 pixels, so that each new pixel covers a linear size of ~ 15.7". This is a good 
compromise between the resolution needed to smooth the very local background variations and at the same time to sample 
the variations due to the different emission from the cluster; 

The exposure maps of the same size and binned by 32 x 32 pixels were constructed at energies of 1.0 keV (soft) and 4.0 keV 
(hard). The values corresponding to each image bin (32x32 original pixels) within the region of the cluster allowed the calculation 
of the local average background per original pixel and the median exposure at the position of the bin. The last step in calculating 
Sum is an estimate of A. To do this it is necessary to study the increasing apparent size of the detected sources as a function of 
the off-axis angle. The size of a detected point source equals, by definition, the size of the Chandra PSF, which depends largely 
on the source's angular distance from the optical axis and, to some extent, on source energy. We studied and interpolated the 
increasing apparent size of the detected sources as a function of the off-axis angle separately in the soft and hard energy bands. 
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Fig. C.l. Sky coverage (area covered vs. flux limit) for the soft (solid line) and hard (dashed line) energy bands computed for the 
area covered by the 18 clusters of the sample. 



For all the clusters we estimated the flux limit corresponding to each region. The sky area available at a given flux limit is then 
simply the sum of all the regions whose minimum detectable flux is lower than S/,,„. The total area covered by the "cluster survey" 
is ~ 0.083 deg 2 . In order to prevent incompleteness effects at the faint end of the source counts, we considered only those sources 
for each cluster with a flux larger than the flux corresponding to 20% of the total sky area covered by each cluster. Figure ICTI 
shows the sky coverage for the eighteen clusters of our sample computed for the two energy bands. The steplike features visible 
in the sky coverage are given by the 20% cut described above. Note that the soft fluxes used to derived the soft-band sky coverage 
are unabsorbed fluxes, corrected for the the Galactic hydrogen column density along the line of sight of each field. 



